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Abstract 


In this research, the application of response surface methodology with central composite design (RSM-CCD) to optimize 
transesterification of palm oil in the presence of zeolite-A as catalyst was investigated. Zeolite-A was synthesized from rice husk 
silica (RHS) and food-grade aluminium foil using hydrothermal method and then characterized using XRD and SEM. The 
synthesized zeolite was then applied to an optimized transesterification reaction using response RSM with three factorial levels, for 
three variables including methanol to oil ratio, catalyst load, and reaction time. The experimental results indicate that the yield of 
99% was achieved at optimum conditions of methanol to oil volume ratio of 6, catalyst load of 9.6%, and reaction time of 4.3 hours. 
The results of experiments and predicted results based on the RSM model are in agreement as shown by the p-value less than 0.05 
at a confidence level of 95%. 
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1. INTRODUCTION biodiesel production. Several examples of vegetable 
oils that have been reported in the literature are 
Awareness of the depletion of petroleum waste cooking oil [1][2], coconut oil [3], Ricinus 


reserves as a traditional source of fossil fuels for 
centuries and the environmental impacts that arise 
from fossil fuel application is a driving factor for 
increasing the use of renewable energy sources. Of 
particular importance are biomass-derived fuels 
(biofuels), biodiesel is one of them that has been 
used as a commercial fuel although still in the form 
of a blend with fossil diesel at a certain proportion. 
Biodiesel is a practical term for the product of 
transesterification, which is a reaction of vegetable 
oil or animal fat with simple alcohol, primarily 
methanol. In accordance with the chemical 
composition of vegetable oil, which is composed of 
free fatty acids (FFAs) and triglycerides, biodiesel 
production involves esterification of FFAs and 
transesterification of triglycerides. However, since 
glycerides are the prominent components of 
vegetable oil, transesterification is the main issue in 
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communis oil [4], rubber seed oil [5][6], and 
soybean oil [7]. Concerning this reaction, biodiesel 
is also known as a mixture of fatty acid methyl 
esters (FAMEs). 

Regardless of its existence which has reached 
commercial status, the price of biodiesel is still 
higher than that of petrochemical diesel, reflecting 
the need for further effort to reduce production 
costs. In this respect, one of the production 
components that has been emphasized is the 
development of heterogeneous catalysts to replace 
homogeneous catalysts that remain as the prime 
catalyst up to the present. In search of effective 
heterogeneous catalysts, various solid catalysts have 
been developed. Several heterogeneous catalysts in 
the form of oxides reported are CaO [2][8], MgO 
[3][9], and calcium methoxide [10][11]. In the form 
of composite, several examples reported in the 
literature are NiO/ZSM-5 [4], MgO-CaO/SiO, [12], 
MgO/Fe203-SiO, [13], MgO/zeolite-A [14], and 
CaO/Si0, [15]. Besides composites, other catalysts 
that have been used are synthetic zeolites, such as 
zeolite-A [16], zeolite-X [17], and zeolite-Y [18]. 

In addition to the type of catalyst, other variables 
acknowledged to influence  transesterification 
reaction are alcohol to oil ratio, reaction time, and 
catalyst load. Considering the role of these reaction 
variables, optimization of the transesterification 
reaction is generally carried out by conducting 
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Table 1. Experimental ranges and levels of the independent variables. 


Variables Symbol Unit 
Catalyst loads A % 
Methanol to oil ratio B - 
Reaction time C hour 


experiments where one variable is varied while the 
others are kept constant. From such experiments, 
the optimum value for the variable being tested will 
be obtained, and this value is then used to determine 
the optimum value for other variables. This method 
has the disadvantage that it does not show an 
interactive effect and the overall effect of the 
variables involved [19]. 
deficiencies, the use of multi-variable statistical 


To overcome these 


techniques is now increasingly explored; one of 
them is known as response surface methodology 
(RSM) [7][20]. RSM is a statistical technique that 
uses multiple regression and correlation analysis as 
a tool to assess the effect of two or more 
independent variables on the dependent variable, 
with the aid of analysis of variance (ANOVA) [6] 
[10][13][21]. This statistical method has been 
applied for the optimization of transesterification of 
different vegetable oils such as papaya oil [22], 
Kusum oil [23], and low-grade cooking oil [24]. 
The present work was conducted to study the 
application of response surface methodology 
(RSM) to optimize the transesterification of palm 
oil in the presence of zeolite-A as a catalyst. 
Zeolite-A was synthesized from rice husk silica 
(RHS) and food-grade aluminium foil using a 
hydrothermal method. Before use, the zeolite was 
using XRD and SEM. The 
then applied to an 


characterized 
synthesized zeolite was 
optimized transesterification reaction using RSM 
with three factorial levels for three variables 
including methanol to oil ratio, catalyst load, and 
reaction time. 


2. MATERIALS AND METHODS 


2.1. Materials and Equipment 
Reagent grade sodium hydroxide (NaOH), nitric 


(5) PANDAWA 


Levels 
-1 0 +1 
5 ie, 10 
3 4 6 
3 4 2 


acid (HNO3), and methanol (CH30H) were 
purchased from Merck. Rice husk and palm oil 
were obtained from local sources in Bandar 
Lampung. The main equipment used was an oven 
(Model Memmert UN Universal 321), furnace 
(Thermolyne Muffle Thermolyne 1100), X-ray 
diffraction instrument (XRD, Xpert MPD), 


scanning electron microscope (SEM, ZEISS EVO 


MA 10), and gas  chromatography-mass 
spectrophotometry unit (GC-MS, QP2010S 
SHIMADZU). 


2.2. Methods 


2.2.1. Extraction of Rice Husk Silica 

Rice husk silica extraction was carried out 
following the previously reported procedure [17]. In 
brief, the process involves several steps including 
mixing 50 g of dry rice husk with 500 mL of 1.5% 
NaOH solution, boiling the mixture for 30 min, 
aging the mixture at room temperature for 24 h, 
filtering the mixture to collect the filtrate (silica 
sol), neutralization of the sol using 10% HNO; 
solution to convert the silica from sol to gel, aging 
of the gel at room temperature for 24 h followed by 
rinsing of the gel with warm distilled water and 
completed by oven drying of the gel for 8 h at 100 
°C, 


2.2.2. Synthesis of Zeolite-A 

Zeolite-A was the 
procedures reported by previous work [25]. To 
commence the experiment, a solvent was prepared 
by dissolving 20 g of NaOH in 250 mL of distilled 
water. A portion of 150 mL of the solvent was used 
to dissolve 30 g of rice husk silica, and the rest to 


dissolve 13.5 g of food-grade aluminum foil. The 


synthesized following 


two solutions were mixed and magnetically stirred 
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Cu-Ka (1.541874 A) 


Z: Zeolite-A (96-900-5742) 
S: Sodalite (96-412-4674) 


Figure 1. The diffractogram of synthesized zeolite-A. 


for 3 h to homogenize the mixture, and then the 
mixture was transferred into an autoclave and 
allowed at room temperature for 24 h (aging 
process). After the completion of this aging process, 
the sample was subjected to crystallization by 
placing the autoclave in an oven at 100 °C for 72 h. 


The solid formed was then removed from the 
autoclave and then washed with distilled water until 
the rinsing water became neutral. The solid was 
oven-dried at 100 °C for 8 h, then calcined at 550 
°C for 6 h, and finally ground into 250 mesh 
powder. 


Table 2. Oil conversion from the experiments and predicted values. 


A: B: 
Rul Catmiystlond (4) Methanol to il ratio 
1 5 3 
2 10 6 
3 75 4 
4 75 3 
5 95 4 
6 5 6 
7 10 3 
8 75 4 
9 10 4 
10 725 4 
W785 4 
122 735 6 
13 10 3 
14 5 3 
Is 5 4 
16 10 4 
75 4 
18 10 6 
ee 3 
20 «5 3 
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Oil conversion (%) 


C: 

Reaction time (h) Actual Predicted 
5 34 31.19 
5 100 98.76 
> 52 55.75 
4 48 53.96 
5 60 55.15 
4 48 47.64 
a) 48 52.02 
3 44 45.99 
3 72 65.11 
5 60 pe fe) 
4 60 53.96 
3 56 54.62 
5 56 52.02 
4 40 37.47 
5 40 38.27 
3 60 65.11 
> 32 38.27 
4 99 101.98 
5 40 44.23 
3 20 21.17 
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Table 3. The results of data analysis using ANOVA. 


Source Sum of Squares 

Model 7130.88 9 

A 3744.58 1 

B 1319.28 1 

C 323.70 1 

AB 188.18 1 

AC 33.46 1 

BC 52.52 1 

A? 23.09 1 

og 0.1520 1 

oy 233.69 1 

Residual 320.07 10 
Lack of Fit 69.40 4 

Pure Error 250.67 6 

Total 7450.95 19 
R’= 0.9570 


2.2.3. Zeolite-A Characterization 

The XRD diffractogram of the synthesized 
zeolite-A was generated by scanning the sample at 
(20) angle in the range of 0—50° on the instrument 
operated with Cu-Ka radiation (A= 1.54 A) with an 
energy of 40 kV and a current of 30 mA. The phase 
of the sample was identified with the aid of Match! 
Version 3.4.2 Build 96 software. The diffractogram 
of the sample was also compared with that of the 
standard zeolite-A provided in the International 
Zeolite Association (IZA) database. The surface 
morphology of the samples was examined using 


Date :12 Apr 2023: UPTLTSIT 
Tine 11:37:19 


Mag= 5.00KX Signal A» SE1 
WO = 11.0mn EHT = 8.00 kV 
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p-value Remarks 


792.32 24.75 < 0.0001 significant 
3744.58 116.99 <0.0001 

1319.28 41.22 < 0.0001 

323.70 10.11 0.0098 

188.18 5.88 0.0358 

33.46 1.05 0.3307 

52.52 1.64 0.2291 

23.09 0.7213 0.4156 

0.1520 0.0047 0.9464 

233.69 7.30 0.0222 

32.01 

17.35 0.4153 0.7928 not significant 
41.78 


Predicted R’= 0.8391 


SEM, with an instrument operated at electron high 
tension (EHT) of 8 kV. 


2.2.4. Transesterification Reaction 

To experiment, an aliquot of 25 mL of oil was 
poured into a 250 mL round bottom flask and then a 
certain volume of methanol, to satisfy the 
predetermined methanol to oil ratio, and a 
predetermined mass of catalyst was added into the 
flask. A small magnetic stirring bar was added to 
the flask. The flask was immersed in a water bath 


paled on top of a hot plate with a magnetic stirrer 


= 10.00KX  SignslA= SE1 
WO = 1.0mm EXT = B.00kV 
7 


Oeste :12 Apr 2023 UPTLTSIT 
Time :114136 


Figure 2. SEM micrograph of synthesized zeolite-A with the magnification of 5,000 (a) and 10,000~ (b). 
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(a) 


Prediction (%) 


Experiment (%) 


Normal % Probability 


Extemally Studentized Residuals 


Figure 3. (a) Plot between experimental results and predicted results, and (b) plot of normal residue. 


and connected with a reflux condenser with water as 
condensing liquid. The reaction was run at 70 °C for 
a specified time. After the completion of the 
reaction, the was cooled to room 
temperature, and the product was filtered into a 
separatory funnel. The funnel was left overnight to 
allow separation between the biodiesel (upper layer) 
and residual oil (lower layer). The volume of the 
biodiesel was measured for the calculation of oil 
conversion using the following equation [5]. 


sample 


% Conversion = 


V, -V, 
- f x 100% 


(1) 


In Equation 1, V; = initial volume of oil and Vs= 
volume of unreacted oil. Identification of the fatty 
acid methyl esters composing the biodiesel was 
conducted using the GC-MS technique, with the aid 
of MS-Library software, NIST 62 and the Wiley 7 
database. 


2.2.5. Experimental Design 

Optimization of the transesterification reaction 
was developed by application of response surface 
methodology with central composite design (RSM- 
CCD) using design expert 13.0 software. In this 


Table 4. Criteria for optimization of reaction. 


Variables Goal 
Catalyst load (%) In range 
Methanol to oil ratio In range 
Reaction time (h) In range 
Biodiesel conversion (%) Maximum 


150 


study, the following second-order polynomial 
model was proposed to predict response (Y) as a 
function of independent variables (Xi) and their 


interactions in Equation (2); 


Y= Bot y BX; + y Bie X7 + yy Bij X:Xj; + € 
i=l i=1 


i=1 j=1 


(2) 


where Y is a response, Bo, Bi, Bi, and Bi, indicate 
the regression coefficients for intercept effects, 
linear effects, quadratic effects, and double 
interactions, respectively. The model reliability of 
the model was evaluated in terms of the correlation 
coefficient, R* obtained from the analysis of 
variance (ANOVA). In this study, Y refers to 
biodiesel yield, and independent variables consist of 
catalyst load, methanol to oil ratio, and reaction 
time. The design of experiment (DOE) with the 
experimental ranges and levels of the independent 
variables is shown in Table 1. 


3. RESULTS AND DISCUSSIONS 
3.1. Zeolite-A Characterization 


Figure | displays the XRD diffractogram and the 
phases composing the zeolite-A identified with the 


Lower limit Upper limit 
5 10 

3 6 

3 ) 

20 100 
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Table 5. The results of model validation at optimum condition. 


A: B: C: 


Catalyst load (%) Methanol to oil ratio Reaction time (h) 


9.6 6 4.3 


aid of Match! Version 3.4.2 Build 96 Software. As 
can be seen in Figure 1, two phases were identified, 
with zeolite-A as the main phase and sodalite as the 
second phase. In addition, the diffractogram also 
displays the existence of an amorphous phase as 
indicated by the presence of peaks with low 
intensity. 

The existence of the sample as a multiphasic 
material is also indicated by the SEM micrograph in 
Figure 2. 

As can be seen, the morphology of the sample is 
characterized by the existence of cubic-shaped 
particles, representing zeolite-A phase, spherical- 
shaped particles coated with rod-like structure, 
representing sodalite, and irregular shape which 
represents amorphous material. 


3.2. Results of Transesterification 

Based on DOE that has been previously 
determined (Table 1), 20 experiments were 
conducted and the percentages of conversion of the 


Conversion (%) 


Conversion (%) 


Prediction Experiment Error 
(%) 2 (%) 
100 99 1 


palm oil obtained from the experiments (actual 
result) were compared to the predicted conversion 
produced with the application of statistical software 
(design expert 13.0 software), as shown in Table 2. 
The data obtained from the experiments were then 
statistically analyzed using the ANOVA test to 
evaluate the effect of variables and the interaction 
between variables applied in the experiments. The 
results are compiled in Table 3. As can be seen in 
Table 3, the p-values of less than 0.05 at the 95% 
confidence level are obtained for the variables A, B, 
C, AB, and C’, suggesting that these variables are 
significant factors in determining the yield of the 
reaction. For the variables AC, BC, A’, and B’ the 
opposite is true. The suitability of the model was 
also evaluated based on several other criteria, 
including insignificant lack of fit, a coefficient of 
determination (R’) that is close to 1, the difference 
between adjusted R’ and predicted R’ is less than 2, 
and adequate precision is more than 4. 

Based on the results of ANOVA as shown in 


A: Amount of 
catalyst (%) 


3 
sew 
SN) 
v 


Figure 4. The 3D plot showing the effect of interaction between catalyst load and methanol to oil ratio (a), 


interaction between catalyst load and reaction time (b), and interaction between methanol to oil ratio and 


reaction time (c). 
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Table 3, the following model equation was 
suggested: 


Y = 66,11 + 20,52A + 11,74B + 6,54C + 6,65AB — 3,45AC + 4,98BC 
— 2,63A2 + 0,4160B?— 11,29C? (3) 


The above model is a second-order polynomial 
quadratic equation, with Y representing biodiesel 
yield, A, B, and C representing catalyst load, 
methanol to oil ratio, and _ reaction time, 
respectively. The positive sign indicates that the 
factor enhances the conversion of the oil into 
biodiesel, while for a negative sign, the opposite is 
true [26]. The suitability of the model is also 
supported by the plot between experimental and 
predicted results as well as the residual normal plot 
as presented in Figure 3. As displayed in Figure 3 
(a), the plot produces a linear line, indicating that 
the predicted values are very close to the 
experimental values. Figure 3 (b) shows that the 
residual points are distributed around the diagonal 
line ideally, indicating that the residuals are 
normally distributed [27]. 

The 3D plot presenting the effect of the 
interaction of the variables applied in 
transesterification (reaction time, ratio of methanol 
to oil, and catalyst load) on biodiesel yield can be 
seen in Figure 4. Overall, the plots in Figure 4 
display the effect of variables and their interaction 
on oil conversion, with the particular purpose of 
obtaining the value of each of the variables which 
leads to the optimum conversion of the oil treated. 

The effect of the interaction between the catalyst 
load and the methanol to oil ratio at a constant time 
of 5 h is shown in Figure 4 (a). The 3D plot shows 
that the increase in the methanol to oil ratio and the 
amount of catalyst is in line with the increase in 
biodiesel yield. Figure 4 (b) shows a 3D plot for the 
interaction effect between the amount of catalyst 


load and reaction time on biodiesel yield, at a 
constant methanol to oil ratio of 6:1. The response 
surface shows that the yield of biodiesel increases 
with increasing loads of catalyst used. Figure 4 (c) 
is a 3D plot for the interaction effect between the 
methanol to oil ratio and reaction time at a constant 
catalyst load of 10%, showing that the yield of 
biodiesel increases with increasing methanol to oil 
ratio. 

In Figure 4 (b) and (c) it can be observed that 
there is a curve associated with the reaction time, 
which indicates that increasing reaction time causes 
an increase in biodiesel conversion to the optimum 
point, followed by a decrease when the optimum 
point was exceeded. In this respect, it can be 
implied that reaction time only slightly influences 
the biodiesel yield. This also means that the 
interaction between the catalyst load and methanol 
to oil ratio is not significant in increasing biodiesel 
yield. This is in accordance with the ANOVA 
results (Table 3) which show that the interaction 
between catalyst load and reaction time as well as 
the interaction between reaction time and methanol 
to oil ratio are not significant in influencing 
biodiesel formation, as indicated by a p-value 
higher than 0.05. The optimization criteria for the 
transesterification reaction applied in this study are 
presented in Table 4. 

Table 5 shows the optimum conditions with 
predicted and validated conversion results, with a 
relatively small error percentage of 1%. This shows 
that the proposed statistical model is suitable for 
optimizing the process of transesterification of palm 
oil with zeolite-A catalyst. 

Many workers have reported the application of 
different zeolite-based catalysts for biodiesel 
production from different feedstocks. Several 
examples are presented in Table 6. As can be seen 
in Table 6, the conversions achieved are varied 


Table 6. Zeolite-based catalysts applied for biodiesel production. 


Catalyst Feedstock Optimum conversion (%) Ref. 


ZnO-B-zeolite Waste cooking oil 


Zeolite-Y Waste cooking oil 
ZSM-S5 Ricinus communis oil 
Zeolite- B Waste sunflower oil 
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84.10 [28] 
84.44 [18] 
96.00 [5] 
66.30 [29] 


32,889,014 
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TIC 


Figure 5. GC-chromatogram of biodiesel produced from transesterification of palm oil. 


depending on the types of catalyst and feedstock 
treated. Compared to the examples presented in 
Table 6 the conversion achieved in this study is 
relatively higher. 

A typical example of GC chromatogram of the 
biodiesel produced is shown in Figure 5. The 
chromatogram is characterized by the presence of 8 
well-separated peaks with varied intensities. This 
feature suggests that there are 8 compounds 
composing the sample. 

The components were then identified with the 
aid of mass spectrometry database (WILEY7.LIB), 
and the results are presented in Table 7. As can be 
seen in Table 7, the eight components of the sample 
are methyl esters that correspond to the fatty acid 
composition of palm oil reported in the literature. 
For example, the composition of palm oil reported 
in the literatures is presented in Table 8. 

As can be seen in Table 8, some composition 
variations between the palm oils from different 
regions are observed. Regardless of such 
variations, it can be seen that palmitic and oleic 


acids are the main components of the oils. 
4. CONCLUSIONS 


In this study, zeolite-A was successfully 
synthesized as confirmed by the results of 
characterization using XRD and SEM. The 
synthesized zeolite-A was found to have high 
activity as a catalyst for the transesterification of 
palm oil into biodiesel, as demonstrated by the 
results of GC-MS analysis which confirms the 
presence of fatty acid methyl esters (FAMEs) as 
components of the biodiesel produced. Application 
response surface methodology with central 
composite design (RSM-CCD) to study the effect of 
three reaction variables (methanol to oil ratios, 
catalyst loads, and reaction times) revealed that the 
variables and their interactions are influencing 
factors for transesterification studied, with the 
statistical model of Y = 66,11 + 20,52A + 11,74B + 
6,54C + 6,65AB - 3,45AC + 4,98BC - 2,63A7+ 
0,4160B* - 11,29C*. Based on the 3D_ plots 


Table 7. Chemical composition of biodiesel produced from transesterification of palm oil. 


Peak no Retention time (min) Componentname Molecularformula Relative amount (“%) 
1 3.461 Methyl caproate C7H,402 0.40 
2 5.488 Methyl caprylate CoH 130 7.76 
3 8.117 Methyl caprate C,,H»O> 6.91 
4 10.219 Methyl laurate C13H2602 41.37 
5 13.405 Methyl myristate C15H3902 20.18 
6 15.405 Methyl palmitate C,7H3402 10.75 
7 17.195 Methyl oleate C19H3602 9.60 
8 17.391 Methyl stearate Cj9H302 3.04 


8) 
Ph) 
22 
3 
= 
> 
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Table 8. Fatty acid components of palm oil reported in the literatures. 


Fatty acid 


Caproic acid 
Caprylic acid 
Capric acid 
Lauric acid 
Myristic acid 
Palmitic acid 
Palmitoleic acid 
Oleic acid 
Stearic acid 
Linolenic acid 
Linoleic acid 


Arachidonate acid 


generated, it was found that the optimum conditions 
for transesterification of palm oil in the presence of 
zeolite-A as catalyst consists of a catalyst load of 
9.6%, methanol to oil ratio of 6:1, and reaction 
time of 4.3 h to achieve the predicted oil 
conversion of 100%. From the experiment, the 
optimum oil conversion of 99% was achieved, 
indicating that the experimental design is in 
agreement with the generated statistical model. 
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